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18.152 Introduction to PDEs, Fall 2011 Professor: Jared Speck 

Class Meeting 7^ 1: Introduction to PDEs 



1. What is a PDE? 

We will be studying functions u = u{x^, x^, - ■ ■ , x") and their partial derivatives. Here 
are standard Cartesian coordinates on M". We sometimes use the alternate notation u{x, y), u{x, y, z), 
etc. We also write e.g. u{r,6,(f)) for spherical coordinates on M^, etc. We sometimes also have a 
"time" coordinate t, in which case t,x^,--- ,x"' denotes standard Cartesian coordinates on M-'^^". 

We also use the alternate notation x^ '= t. 

We use lots of different notation for partial derivatives: 

d 

(1.0.1a) ^^"^ ~ ''^^^ ~ ^^'^^ ^ ^ i ^ n, 

(JX 

(1.0.1b) r. ■ r. ■ = -7r~^ — ^"^ = Ux^xi = didjU, 1 < i,j < n. 

ox^ox^ ax* ox^ 

If i = j, then we sometimes abbreviate didjU '= dfu. If -u is a function of (x, y), then we also write 
Ut — T{ ti. etc. 

Definition 1.0.1. A PDE in a single unknown u is an equation involving u and its partial deriva- 
tives. All such equations can be written as 

(1.0.2) F(m,m^i,--- ,u^.n,M^i^i, ■ ■ ■ ,M^,i...^»jv,a;\x^- ■ ■ ,x") = 0, ii, ■ ■ ■ , e {1, 2, ■ ■ ■ ,n} 
for some function F. 

Here N is called the order of the PDE. N is the maximum number of derivatives appearing in 
the equation. 

Example 1.0.1. u = u{t,x) 

(1.0.3) -d^u + {l + cosu)dlu = 

is a third-order nonlinear PDE. 
Example 1.0.2. u = u{t,x) 

(1.0.4) -d^u + 2dlu + u = t 

is a second-order linear PDE. 



We say that (1.0.4) is a constant coefficient linear PDE because u and its derivatives appear 
linearly (i.e. first power only) and are multiplied only by constants. 
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Example 1.0.3. u = u{t,x) 



(1.0.5) 

is a third-order linear PDE. 



dtu + 2(1 + x^)dlu + u = t 



We say that (1.0.5) is a variable coefficient linear PDE because u and its derivatives appear 
linearly (i.e. first power only) and are multiplied only by functions of the coordinates {t,x). 

Example 1.0.4. u = u{t, x),v = v{t, x) 

(1.0.6a) dtu + 2xdxV = sin(a;^), 

(1.0.6b) dtv - x^d^u - 

is a system of PDEs in the unknowns u,v. 



2. The Goals of PDE (and of this course) 

Suppose that we are interested in some physical system. A very fundamental question is: 

• Which PDEs are good models for the system? 

A major goal of modeling is to answer this question. There is no general recipe for answering it! 
In practice, good models are often the end result of confrontations between experimental data and 
theory. In this course, we will discuss some important physical systems and the PDEs that are 
commonly used to model them. 

Now let's assume that we have a PDE that we believe is a good model for our system of interest. 
Then most of the time, the primary goals of PDE are to answer questions such as the following: 

(1) Does the PDE have any solutions? (Some PDEs have NO SOLUTIONS whatsoever!!) 

(2) What kind of "data" do we need to specify in order to solve the PDE? 

(3) Are the solutions corresponding to the given data unique? 

(4) What are the basic qualitative properties of the solution? 

(5) Does the solution contain singularities? If so, what is their nature? 

(6) What happens if we slightly vary the data? Does the solution then also vary only slightly? 

(7) What kinds of quantitative estimates can be derived for the solutions? 

(8) How can we define the size (i.e., "the norm") of a solution in way that is useful for the 
problem at hand? 



3. Physical Examples 

It is difficult to exaggerate how prevalent PDEs are. We will discuss some important physically 
motivated examples throughout this course. Here is a first look. 

• —d'^u + d^u = wave equation, second-order, linear, homogeneous 

• —dtU + d'^u = heat equation, second-order, linear, homogeneous 

• d'^u + dyU + d^u = Laplace 's equation, second-order, linear, homogeneous 

• d'^u + dyU + d'^u = f{x,y, z) Poisson's equation with source function /, second-order, linear, 
inhomogeneous (unless / = 0) 

• idtu + diu = Schrddinger's equation, second-order, linear, homogeneous 

• Ut + Ux = 0, transport equation, first-order, linear, homogeneous 

• Ut + uUx = 0, Burger's equation, first-order, nonlinear, homogeneous 
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E = {Ei{x,y,z),E2{x,y,z),Es{x,y,z)), B = {Bi{x,y, z), B2{x,y, z), B3{x,y, z)) are vectors in 

(3.0.7a) 9tE - V X B = 0, V • E = 0, 

(3.0.7b) 9tB + V X E = 0, V • B = 

"Maxwell's equations" in a vacuum (i.e., matter-free spacetime), first-order, linear, homogeneous. 

4. Linear PDEs 

Before we dive into a specific model, let's discuss a distinguished class of PDEs that are relatively 
easy to study. The PDEs of interest are called linear PDEs. Most of this course will concern linear 
PDEs. 

Definition 4.0.2. A linear differential operator £ is a differential operator such that 

(4.0.8) £{au + bv) = a£u + hCv 

for all constants a, 6 e R and all functions v. 

Renicirk 4.0.1. The notation was introduced out of convenience and laziness. The definition is 
closely connected to the superposition principle. 

Example 4.0.5. C =^ —d"^ + if' — x'^)dl is a linear operator: Cu = —d^u + (t^ — x^)dlu 

Example 4.0.6. u = u{x,y), Cu = d^u + u^dyU does NOT define a finear operator: 
C{u + v) = dl{u + v) -\- {u + vYdy{u + v) ^ d^u + v?dyU + d^v + v'^dyV = Cu + Cv 

Definition 4.0.3. A PDE is linear if it can be written as 

(4.0.9) £m = /(,t\--- ,x") 

for some linear operator C and some function / of the coordinates. 

Definition 4.0.4. If / = 0, then we say that the PDE is homogeneous. Otherwise, we say that it 
is inhomogeneous. 

Example 4.0.7. u — u{t,x) 

(4.0.10) dtu - (1 + cos t)dlu = tx 

is a linear PDE. 

Here is an incredibly useful property of linear PDEs. 
Proposition 4.0.1 (Superposition principle). ■ ■ ■ ,mm are solutions to the linear PDE 

(4.0.11) £w = 0, 

and ci, • • • , Cm G R, then YldLi '^i'^i ^■^ (^^^o a solution. 
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Proof. 

(4.0.12) £ aui = ^ Ci Cui = 0. 

1=1 i=l 

□ 

Remark 4.0.2. This shows that the set of all solutions to Cu = is a vector space when C is 
linear. 

As we will see in the next proposition, inhomogeneous and homogeneous linear PDEs are closely 
related. 

Proposition 4.0.2 (Relationship between the inhomogeneous and homogeneous linear 
PDE solutions). Let Sh be the set of all solutions to the homogeneous linear PDE 

(4.0.13) Cu = 0, 

and let uj be a "fixed" solution to the inhomogeneous linear PDE 

(4.0.14) Cu = f{x\--- ,x''). 



Then the set Sj of all solutions to (4.0.14) is the translation of Sh by uj : Sj = {uj + uh \ uh ^ 
Sh}. 



Proof. Assume that Cui = /, and let w be any other solution to (4.0.14), i.e., Cw = f. Then 



C{w — uj) = / — / = 0, so that w — uj G Sh- Thus, w = uj + {w — ui) , and so w ^ Si 

^ V ' 

belongs to Sh 

by definition. On the other hand, if w G S"/, then w = uj + uh for some uh G Sh- Therefore, 



Cw = C{ui + Uh) = Cuj + Cuh = f + ^ = f- Thus, w is a solution to (4.0.14). □ 

5. How TO SOLVE PDEs 

• There is no general recipe that works for all PDEs! We will develop some tools that 
will enable us to analyze some important classes of PDEs. 

• Usually, we don't have explicit formulas for the solutions to the PDEs we are 
interested in! Instead, we are forced to understand and estimate the solutions without 
having explicit formulas. 

The two things that you typically need to study a PDE: 

• You need to know the PDE. 

• You need some "data." 

6. Some simple PDEs that we can easily solve 

6.1. Constant coefficient transport equations. Consider the first-order linear transport equa- 
tion 

(6.1.1) adxu{x, y) + bdyu{x, y) = 0, 

where a, 6 G M. Let's try to solve this PDE by reasoning geometrically. Geometrically, this equation 
says that Vu ■v = 0, where Vu '= {d^u, dyu) and v is the vector (a, b) G M^. Thus, the derivative of 
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u in the direction (a, b) is 0, which imphes that u is constant along hnes pointing in the direction 
of (a, 6). The slope of such a line is ^. Therefore, every such line can be described as the set of 
solutions to bx — ay = c, where c G M. Since u is constant along these lines, we know that u is a 
"function that depends only on the line c." Therefore u{x, y) = f{c) = f{bx — ay) for some function 
/• 

In order to provide more details about u, we would need to prescribe some "data." For example. 



if it is known that u{x, 0) 



X 



then x^ = f{bx). Thus, /(c) = b ^c^, and u{x,y) = {x — b ^ay) 



In the future, we will discuss the kinds of data that can be specified in more detail. As we will see, 
the type of data will depend on the type of PDE. 

6.2. Solving a variable coefRcient transport equations. With only a bit of additional effort. 



the procedure from Section 6.1 can be extended to cover the case where the coefficients are pre- 
specified functions of x,y. Let's consider the following example: 



(6.2.1) 

Let P denote a point P 



ydxU + xdyU = 0. 



{x,y), and let V denote the vector V = {y,x). Using vector calculus 
notation, (6.2.1) can be written as Vu{P) • = 0, i.e., the derivative of m at P in the direction 



of V is 0. Thus, equation (6.2.1) implies that u is constant along the curve C passing through P 



that points in the same direction as V. This vector can be viewed as a line segment with slope ^. 
Therefore, if the curve C is parameterized by x — ?■ {x,y{x)) (where we are viewing y as a function 
of x along C) then C has slope ^, and y is therefore a solution to the following ODE: 



(6.2.2) 



dx 



X 

y' 



We can use the following steps to integrate (6.2.2), which you might have learned in an ODE 
class: 



(6.2.3) 
(6.2.4) 



(6.2.2) 



y 



dy_ 
dx 

2 



r _ 3^ 

2 ~ 2 

Thus, the curve C is a hyperbola of the form {y^ — x"^ 



1 d 



iy' 



+ c, c 



2dx 

constant. 



X 



c}. These curves are called characteristics. 



We conclude that u is constant along the hyperbolas {y'^ 
f{x^ — y"^) for some function /(c). 



X 



c}, which implies that u{x,y) 



We can carry out the same procedure for a PDE of the form 



(6.2.5) a{x,y)dxU + b{x,y)dyU = 0, 

as long as we can figure out how to integrate the ODE 



(6.2.6) 



dy 
dx 



b{x,y) 



a{x,y) 

7. Some basic analytical notions and tools 



We now discuss a few ideas from analysis that will appear repeatedly throughout the course. 
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7.1. Norms. In PDE, there are many different ways to measure the "size" of a function /. These 
measures are called norms. Here is a simple, but useful norm that will appear throughout this 
course. 

Definition 7.1.1 {C'^ norms). Let / be a function defined on a domain C M. Then for any 
integer A; > 0, we define the norm of / on i7 by 

k 

a=0 

where f'^°'\x) is the a*'^ order derivative of f{x). We often omit the symbol Vt when = M. 
Example 7.1.1. 

(7.1.2) ||sin(x)||cr(R) =8. 

The same notation is used in the case that Vt C M", but in this case, we now sum over all 
partial derivatives of order < k. For example, if f2 C M^, then ||/||c2(f7) == sup(^ t/)| + 
suP(^,j,)6Q l'9x/(a;,?/)|+sup(^^^)gf^ \dyf{x,y)\+s\x^(^^^y^^^ \dlf{x,y)\+sn^^^^^y)^^ ?/)|+sup(^^y)gf^ \dlf\x,y)\. 

If / is a function of more than one variable, then we sometimes want to extract different informa- 
tion about / in one variable compared to another. For example, if / = /(t, x), then we use notation 
such as 

1 2 

(7.1.3) 11/11^.,.^=^^^ sup \dtf{t,x)\ + Y, sup \d:f{t,x)\. 

Above, the "1" in C^'^ refers to the t coordinate, while the "2" refers to the x coordinate. 

The next definition provides a very important example of another class of norms that are prevalent 
in PDE theory. 

Definition 7.1.2 (L^ norms). Let 1 < p < oo be a number, and let / be a function defined on a 
domain C M". We define the norm of / by 

(7.1.4) ll/IUnn) = (^|/(x)rrf"x 

We often write just instead of Lp(M"). 

II ■ I|lp(c) has all the properties of a norm: 

• Non-negativity: ||/||lp(q) > 0, ||/||lp(q) = <^=^ f{x) = almost everywher(]^ 

• Scaling: ||A/||iP(n) = |A| ||/||iP(Q) 

• Triangle inequality: \\f + g\\LP{n) < \\f\\LP{n) + Wghpin) 
Similarly, || . ||cfe(n) also has all the properties of a norm. All of these properties are very easy to 
show except for the last one in the case of || ■ ||LP(n). You will study the very important case p = 2 
in detail in your homework. 




'Almost everywhere" is a term that would be precisely defined in a course on measure theory. 
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7.2. The divergence theorem. A lot of PDE results are derived using integration by parts (some- 
times very fancy versions of it), which provides us with integral identities. This will become more 
apparent as the course progresses. Let's recall a very important version of integration by parts from 
vector calculus: the divergence theorem. We first need to recall the notion of a vectorfield on M". 

Definition 7.2.1 (Vectorfield). Recall that a vectorfield F on C M" is an M"— valued (i.e. 
vector-valued) function defined on Q. That is, 

(7.2.1) F : ^] ^ M", 

F(x\--- ,x") = (Fi(x\--- ,x-),--- ,F-ix\--- ,x-)), 

where each of the are scalar- valued functions on M". 

We also need to recall the definition of the divergence operator, which is a differential operator 
that acts on vectorfields. 

Definition 7.2.2 (Divergence). Recall that V-F, the divergence of F, is the scalar-valued function 
on M" defined by 

n 

dcf 



(7.2.2) V ■F = J2diF\ 



i=l 



We are now ready to recall the divergence theorem. 

Theorem 7.1 (Divergence Theorem). Let Q cM.^ be a domaii^with a boundary that we denote 
by dQ. Then the following formula holds: 

(7.2.3) I V ■¥{x,y,z)dxdydz = I F(a) ■ N((t) da. 

Jn JdQ 

Above, N(cr) is the unit outward normal vector to dQ, and da is the surface measure induced 
on dQ. Recall that if dQ C M.^ can locally be described as the graph of a function (f){x,y) (e.g., 
dQ = {{x,y,z) \z = (j){x,y)}), then 



(7.2.4) rfa = V 1 + I V(/)(x, y)\'^dxdy, 

where V0 = {dx(f),dy(f)) is the gradient of (p, and |V0| \/ {dx<\>Y' + {^yCpY is the Euclidean length 
ofV<p. 

Remark 7.2.1. The divergence theorem holds in all dimensions, not just 3. In dimension 1, the 
divergence theorem is 

(7.2.5) / —F{x)dx = F{b)-F{a), 

J[a,b] dx 

which is just the Fundamental Theorem of Calculus. 



Throughout this course, a domain is defined to be an open, connected subset of ] 



MIT OpenCourseWare 
http://ocw.mit.edu 



18.152 Introduction to Partial Differential Equations. 
Fall 2011 



For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms . 



